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1 Introduction 



The quasi-one-dimensional spin ladder models have attracted much attention in recent years; 
besides the mathematical interest, these models are believed to describe many real magnetic 
materials invented relatively recently 

In connection with this the construction and exact solution of exactly integrable ladder 
models became very actual. In particular the Bethe Ansatz exact solution can provide a deep 
insight into the properties of these models. Recently in a series of articles [HISISISI, it have 
been proposed a technique for constructing integrable models where the model parameter A (the 
anisotropy parameter of the XXZ and anisotropic t — J models) have a staggered disposition 
of the sign along both chain and time directions. Due to the staggered shift of the spectral 
parameter also considered there, these models have led to Hamiltonians formulated on two leg 
zig-zag quasi-one dimensional chains. 

The technique is based on an appropriate modification of the Yang-Baxter equations (YBE) 
( jni 13 IHl ) for the i?-matrices which are the conditions for commutativity of the transfer matrices 
at different values of the spectral parameter. The transfer matrix is defined as the staggered 
product of i?-matrices, which have staggered sign of the anisotropy parameter A and staggered 
shift of the spectral parameter of the model along the chain. The shift of spectral parameter in a 
product of i?-matrices was considered in earlier works of several authors [HI El , by use of which 
they have analyzed the ordinary XXZ model in the infinite limit of the spectral parameter. 
Later, in the article ^,11^, this modified transfer matrices has been used in order to construct 
integrable models on the ladder. 

Various integrable models on ladders were also constructed in articles ^1 El El E] • 
The above described construction essentially differs from those due to the inhomogeneity of 
anisotropy parameter along the chain. 

The subject of present paper is to study the exactly soluble spin-1/2 two-leg zig-zag ladder 
model proposed in 0, in the thermodynamic limit. We use the technique of Thermodynamic 
Bethe Ansatz developed by M. Gaudin [221 and M.Takahashi and M.Suzuki j26j for the Heisen- 
berg XXZ model. This method is based on the study of Bethe Ansatz equations (BAE). 
Another powerful method in theoretical study of spin ladder models which we partially use is 
the bosonization method |171 1181 ITU] . The bosonization have been used very often to find out 
and to understand many physical aspects of these systems . This method yields a low-energy 
effective field theory of a spin system, which can be treated by powerful methods of quantum 
field theory. 

In the model under consideration, the i?-matrix is a product of two samples of XXZ R- 
matrices with the staggered sign of anisotropy parameter A and the shift of spectral parameter 
by a new model parameter. The Hamiltonian 
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describes the interaction of 2N spins located on the sites of zig-zag ladder (see Fig^, in the 
external "magnetic" field H along the z-axis. Parameterizations of coupling constants Ji and 
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Figure 1: Zig-zag ladder chain 
J2 are given by the two combinations 

1. A ~ real, 6 ^ imaginary, then < | Ji| < | J2I < cxd. 



2. A ~ imaginary, ~ real, then < | J2I < | Ji| < 00. 



(2) 



which fix the Hamiltonian to be Hermitian and give all the possible (imaginary) values of Ji 
and J2. 

By construction, this model can be solved exactly through the Bethe Ansatz technique 
Starting from the ferromagnetic eigenstate with all spins up, each M- particle eigenstate is ob- 
tained by adding M magnons which are parameterized by unequal complex rapidities xi, ..,xm, 
the roots of the Bethe Ansatz equations [2j 
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The corresponding state has total spin projection Sz = N — M and energy 
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This type of BAE arise also in the study of lattice limits of one-dimensional field theoretical 
models j2Ul 1211 ^FR I23j . as well as in some exactly solvable zig-zag ladder models considered in 
jlll I24j . But the similarity is restricted only by the Bethe equations, because for all of these 
models the bare energy h{x) is even function of the rapidity x in contrast of our case, where h 
is apparently odd function of x. As we shall see, this makes the physical pictures very different. 

The paper is organized as follows. In section 2 the BAE and the classification of their 
solutions in different sectors of the model are given, in the thermodynamic limit N ^ 00. The 
spectral equations provided that the system is in the statistical equilibrium state are derived. 
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In section 3, the spectral equations are considered in zero "magnetic" field H = 0. We have 
reduced the system of spectral equations in zero temperature limit T ^ to a single integral 
equation and sketched out the structure of ground states and lowest lying excitations of the 
model, for two different sectors and We shown that the spectrum is gapless in and 
has both gapped and gapless excitations in @. Two special limits Ji = 0, J2 7^ and Ji = J2 
are considered in section 4. We have shown that when Ji = 0, in the thermodynamic limit our 
model is equivalent to two species of Heisenberg XX Z models with different signs of hopping 
constants. We studied the interaction between the two species, which is of next to leading order, 
by use of both the finite size correction and the bosonization methods. We recognized that the 
limit Ji = J2 inherits some features from both |Ji| < IJ2I and |Ji| > IJ2I regions. We have 
summarized our results in the last section. 

2 The Bethe Ansatz equations in the thermodynamic hmit and 
the spectral equations 

In what follows we shall be interested in the thermodynamic limit of the theory. Namely, the 
size of the system N and the number of elementary excitations in the energy interval [e, e + de], 
M(e), will go to infinity by the finite ratio 

N ^ 00, M(e) ^ 00, M(e)/A^ = const 

Then, according to [^511101) general solutions of © are grouped in so-called strings. The char- 
acter of the string solution strongly depends on whether the proper value of complex parameter 
A is real or imaginary. Thus in two regions and ^ the model exhibits different properties 
which we shall analyze separately. 

In spite of possible shortcomings of the string hypothesis ( see for instance |27| ), we base 
our analysis on the string solution approach, which we believe renders correct results in this 
thermodynamic regime, as it is in the case of Heisenberg spin chains. 

It is worth to mention that we take the parameter 9 to be finite. In the absence of magnetic 
field, all the spectral functions typically go to zero at (finite or infinite) Fermi points. Then 
we expect that in the region ^ there will be no effective restrictions on the variation range of 
spectral parameters, and the spectrum will be gapless. 

In this section we shall bring the basic properties of Bethe Ansatz exact solution. This 
solution will allow to derive the spectral equations of the theory which describe the spectrum of 
physical excitations and provide a systematic method for constructing the physical vacuum and 
the equilibrium properties of the system. 

2.1 The case of real A and imaginary 9 (|Ji| < | J2I) 

2.1.1 The properties of Bethe Ansatz equations 

In this case we change the parameter 6 ^ iO in order to deal with real 9. 

Equations © and eigenvalues @ have periodic property with respect to the shift of param- 
eters Xi —>■ Xi + 2pQi, where the parameter = f is introduced. In the thermodynamic limit, 
string solutions have the form |26j 

xl'X = x]^^±+i{n + l-2k)+i^{l-v), (mod2po«), A;=l,2,..,n; (5) 
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characterized by the parity v = ±1, the real part ( the center of string ) ^ and the order n. 
The allowed values of n are restricted by the conditions 



V sin[/cA] sin[(n — k)X\ > 0, 



(6) 



which follows from the normalizability conditions for the wave functions j28j . 

Due to the symmetry A — > vr — A of the eqs. we shall first consider the case A < 7r/2. 
Solutions of © for arbitrary values of parameter pQ are given in The starting point is the 
representation of po in terms of continued fraction 



Po = ^'1 + 



V2 + 



1 ' 
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where Vi are integers. Now, in order to describe strings which satisfy one defines series of 
real numbers pi and series of integers rrii and yi as follows: 



po = 7r/A, pi = 1, Pi = pi-2 - Pi-iVi-i, 

i 

mo = 0, mi = ^i/fc, 

k=l 

y-i = 0, yo = 1, yi = vi, yi = yi-2 + i^iyi-i- 
The order and parity of all strings can be expressed by these series: 



(8) 
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Vi-i + (j - m)yi for mi<j< TTT-i+i, j = 1, 2, 
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-1, Vj = exp{Tri[{nj - l)/po]), 



where " [ ]" is the Gauss' symbol. Suppose that there are Mj bound states of parity Vj and order 
n,-. Then we write the coordinate of the center Eqs. can be written in terms of 

the centers of strings 

OO Mn 
71=1/3=1 

Here the function Pj{x) is the bare j— string momentum 



X 



(10) 



Pj{x) = fix + 9/X,nj,Vj) + fix - e/X,nj,Vj), 
and the bare scattering phase between the j— string and the n— string equals 

Qjkix) = '^[fix,p,VjVk) + fix,p + 2,VjVk)], 



(11) 



(12) 



where the sum is over p = |?7-j — n^l, |nj — n^l +2, .., nj +n/c — 2. In eqs. (|11() and H12|) the function 
sinh(|(x-i/c)-if(l-f))^ 



fix,k,v) = 
is introduced. 
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The total momentum of a given configuration is P = ^Pkix). The energy can be expressed 

k 

as E = A'^ hk{x), where 
k 

hk{x) = ^[f{x - 0/X,nk,Vk) - f{x + e/X,nk,Vk)] (14) 

is the bare-particle energy, and the constant 

cosh(20) — cos(2A) tan A 



A 



sinh^ e A 



is introduced. 

The quantities Jq are integers or half-integers then N is odd or even respectively. The set 
{Ja} defines a unique solution {xa} of A set of numbers Ji omitted in {Ji} defines the 
holes in the distribution of Xa. In the thermodynamic limit, the string "particles" and "holes" 
are continuously distributed with densities pj and pj. Then eqs. © in the leading order in 1/A'' 
acquire the form of integral equations as follows: 

b,{x) = + pj{x)) + Y,Tjk * Pk{x), (15) 

k 

where 

and ^ ^ 

bk{x) = ^^[/(^ + ^/Knk,Vk) + fix - e/X,nk,Vk)]. 
The convolution "*" is defined as 

+00 

/ * 9{x) = / /(x - y)g{y)dy. 



The sign-factor (— 1)^'0) = sign(Ja — Jl^_i) indicates the signature of the derivative f'{x,j,Vj] 
It equivalently can be defined from the inequality (see Ref. 

The energy and momentum are of the form 



Following to [2S^, we can express the entropy by 



00 °? 

E/N = Yj {-^hjix) + 2njp,oH) pj{x)dx - poH, (16) 

•J=l-oo 
00 

00 „ 

P/N = Y1 J Pj{x)Pj{x)dx. (17) 



00 °? 

^/^ ^12 ^^Pi + ^"^^Pi + Pi^ ~ Pj Pi ~ Pj Pj^ 
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To treat the case A > 7r/2, one can use the above mentioned symmetry A ^ vr — A and replace 
A by vr — A in H1U|) - ()15|) . This replacement preserves the previous classification of strings, with 
the new parameter po = ^p^. Eqs. retain their form, but the bare energy and momentum 
change their signs in and H14|) : 

(x. A) — > — vr — A), A) — > — /ij(x, vr — A). (19) 



2.1.2 The spectral equations 

Let the system is in the state characterized by densities Pj{x) and Pj{x). The equilibrium 
dynamics of the system at some temperature T can be extracted by minimizing the free energy 
F = E — TS with respect to independent pj. This yields the following non-linear integral 
equations for r/j = Pj/pj (the main spectral equations): 

^ oo 

Inr?,- = - {Ah, + 2n,^i^H) + ^(-1)^'=)^,. , * in(l + r?, i). (20) 

k=l 

It is more convenient to deal with a symmetric integral operator 

A,k{x -y) = (-1)^(^)<5,, fc<5(x -y) + T,- fc(x - y). 
Its inverse operator has the form 

= (-1)'^^ -y)- Si[{l - 26m,^,j)6j.i^i + 5j+i^i]) , mi_i < j < rm - 2, 

(21) 

Aji^ = (-1)*'*'^ -y) - Si{l - 26m,_i,j)Sj-i,i - didj^i - Si+i6j+i^i) , j = rrii - 1, 



where 



00 



1 1 J f duo e ^"^^ cosh.[{pi — pi^iju:] 

' nx 



4pjCOsh(^)' ^«(^) y 27r 2 cosh(pjCi;) cosh(pj+ia;) 



(22) 



This operator acts on hj and hj in very simple way: 

^^7/ * hi = 27r5j- 1 (si(x - e/X) - si{x + e/\)) = eo{x). 

J] ^7/ * h = 6j,i (si(x + 9/X) + si{x- e/X)) , 
I 

The free energy can be expressed in transparent form 

00 

F = -T j {si{x + e/X) + si{x-e/X))ln{l + r]i)dx. (23) 

—00 

Here we used (|lf)j) . (PH]) and the properties (|221)- 
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In the zero temperature limit T ^ the system goes to its vacuum state. For further analysis 



of this limit it is useful to introduce the so-called pseudo-energies by rjj{x) = exp 
transform 1)20(1 in terms of these functions as follows: 

oo 

Ahj{x) + 2njfioH - fc * e^{x), 



T 



k=l 



or, using the properties of inverse operator (|22() . 

ei(x) = Ae(j{x) + si* 62 (x), 



Ejix) = Si* [(1 - 26m,_-,,j)ej_^{x) + e++i(x)], mj_i < j < - 2 



ej{x) = Si * (1 - 2(5m^_j +di* e+(x) + Sj+i * e/+i(x), j = - 1, 



j^oo rij 

Here the dagger means positive parts of corresponding functions: 

ej{x), if £j{x) > 0, 

^ 0, if Sjix) < 0, 
Prom eqs. (|16|) and ()24|) one can represent the energy in a transparent form 



E — En + E, 



~r -^exci 



where 



E^/N =^(-l)'^(^) j ejix)bjix)dx 

[si{x + 9/X) + si{x- e/X)] ef{x)dx, 



fix. 



and 



(24) 



(25) 



(26) 

(27) 
(28) 



By definition, (|28() is non-negative. Therefore, i?mm = -E'o is the energy of vacuum configuration, 
which can be realized then the following two conditions are satisfied: 



eUx)pj{x) = 0, e- {x)pj{x) = 0. 



(29) 



The expression for vacuum energy (|27() follows from H15|). (|21() . (|24l) and agrees with (|23() in the 
limit r ^ 0. 



2.2 The case of imaginary A and real 6 (| Ji| > | J2I) 
2.2.1 The properties of the Bethe Ansatz equations 

In order to ease notations and deal with real spectral parameter and real A, we change x — > j(p 
and A — > iX. The corresponding Bethe Ansatz equations and all the physical quantities are 
periodic in each 99^ with a real period 2tt. 
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The string solutions are of the form 



^''^ - ipl+i\{n + l-2k), (mod27r), A; = l,2,..,n; 



(30) 



where (p"^ is the real part ( the center of string ) and the order n is not bounded from above |25j . 
The bare particle momentum Pni'-p) and energy are defined like in and (|14() : 



Pn{x) = f{(p + e,n) + f{(p - e,n), 
hk{x) = -^[f{v + e,n)-f{^-9,n)]. 



(31) 
(32) 



Here the function / is defined as 



-i In 



sin ^(v — i^k) 
sin ^{^ + i^k) 



(33) 



Due to the mentioned periodicity, one can restrict the interval (pk £ ["'^i "^1 once fixing the 

logarithmic branch in (|33j) . 

The bare two-particle scattering phase Qjk is also given by (fT^ with the function (1321) • 
The densities of the n-strings and holes Pni'p) and Pni'p), defined in the region [— 7r,7r], obey 

the following integral equations: 



where 



(34) 



The symbol "*" means 



— TT 

The energy, momentum and entropy of a given configuration Pnip), Pni'-p) are 

oo ^ 

E/N = ^ {AK{^) + 2npoH) Pn{^)dp - poH, 

oo J 

P/N = Y1 / Pn{p)Pn{^)dp. 
n=l 

— TV 
OO J 

S/N = ^ / {{pn + Pn) ln(/9„ + Pn) " Pn In Pn " Pn In Pn) dip. 



(35) 
(36) 
(37) 
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2.2.2 The spectral equations 

In this case, the same steps as in previous one are appropriate. First we minimize the free energy 
and get the following spectral equations: 

^ oo 

lnr?„ = - {AK + 2nfioH) + ^ T„,^ * ln(l + r]-^). (38) 

m=l 

Then we pass to zero temperature limit, where ()38() is equivalent to the following integral 
equations for pseudo-energies rjn{^p) = exp ^ ^"jf^ ^ : 

oo 

£+((/?) = Ahn{<^) + 2n/ioi? - ^ An^rn * ^mi'^)^ (39) 

m=l 

where we introduced the operator 
Its inverse will have a form 

Anlrni'P) = ^n,m5{<f) - ^D{(f){6n,m-l + 5n,m+l), (40) 



with 



2cosh(A;A) 

k 

which is positive valued function with period 27r. Inverting the system ()39() . one gets: 

ei{ip)=Aeo{ip) + ^D*et(.^), (41) 
£„((/.) = * + 4+1(95)) , n > 2 

hm ^ = 2/xoi/. 

n— >oo fi 

Here the property 

An!m * hM = <5i,neo(^) = h,n {D{^ + 9)- D{ip - 9)) 

have been used. 

This system of coupled equations is simpler than the previous system (|25|) . In particular, 
one can immediately conclude that the functions e„ are non-negative for n > 2. This means 
that the vacuum is formed by 1- strings only. 

For the energy of given configuration p„, p„,, the equations analogous to ()26() - ()28() are valid: 

E = Eq + Eexci 
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where 



Eo/N = Y,I e-i^)bni^)d^ (42) 

n 

TT 

= -^l i^i^ + ^) + - ^)] (43) 

— TT 

^exc/iV = / [4(</')Pn(</') - e-{ip)pniv)] dip. (44) 

3 Ground state and excitation spectrum in zero "magnetic" 
field 

Here we shall study the structure of the ground state and excitation spectrum at H = 0. 
First consider the case of real A, in the region < A < 7r/2. By definition, 



Vjix) = pj{x)/pj{x) = exp 



T 



In the limit T — > survive only the states whose rapidities obey £j{x) < 0. Relations ()29() also 
indicate that, in the ground state, the occupied rapidities obey ej < 0. Let us denote this region 
by Aj (the Dirac sea): 

£j{x) < 0, X G Aj. 

According to this, the structure of the ground state is governed by the system of coupled equa- 
tions (123). 

It can be shown that in zero "magnetic" field H = eqs. H25|) have solutions with the 
following properties: 



j G,ix-y)eiiy)dy>0, 2 < j < mi - 1, (45) 



ei>0 
oo 



' 7 27r cosh(po-2)a; " -J- ^ 

— oo 

em^x) = -S2 * £^1-1(2;) < 0, (46) 

ej(x) = 0, i>mi. (47) 

These relations constitute the solution of H25|) in terms of unknown function £i{x). Inserting 
the integral representation (|45() for non- negative £2 into the first equation of H25() . we arrive to 
the following integral equation with respect to ei(x): 

£i{x) = A£q{x) + j R{x-y)£i{y)dy, (48) 

£1>0 

00 

R{x) = s,*G2{x)= I ^e— cosh[(po - 3).;] 

^ ' ^ ' J 2-K 2cosh[(po-2)a;]cosh[a;] - 
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Here one of course needs to solve eq. . One could assume that the Wiener-Hopf method 
would be effective. This method would nicely work in the case when the region ei > is 
spanned from infinity to some finite point. Our analysis shows that this is not the case. Then 
We employed both perturbative Wiener-Hopf and numerical methods. In order to be shorter, 
we will present the detailed analysis and related results elsewhere. Meanwhile we bring here 
some speculative analysis which make the general picture rather clear. 

As defined by (|22]) . the so-called disturbance function eq is an odd function with positive 
values at positive arguments, £o{x) > 0, x > 0. This means that ei is also positive at non- 
negative values of x: ei{x) > 0, x > > 0). Further analysis of eq. (|48j) shows that £i{x) 
takes negative values in some region in (— oo, 0). This region which we denoted by Ai, is the Dirac 
sea of 1- strings. Furthermore, Emi is non-positive within the whole real axis — oo < x < -|-oo. 
So the Dirac sea of mi- strings A^i is (— oo,-|-oo). The rest pseudo-energies are non negative. 
Thus the vacuum is filled by the 1 and nii strings: 

Ai C (-00,0), 

A™, = (-oo,+oo), (49) 
A, = 0, 

A = Aie A^^. 

where p^, are vacuum distribution functions, A is the whole Dirac sea. Integrating (|15() for 
j = mi we get that the vacuum configuration is + = N, Mj =0, j 7^ 1, mi, where M^^ 
are numbers of vacuum i-strings. So we have the half-filled band and Sz = 0, i.e., we deal with 
an antiferromagnetic case. 

Let us turn to considerations of excitation spectrum. The possible excited states can be 
formed by adding 1^ excited /c-strings with rapidities z^, a = 1, ..,lk and subtracting Vj vacuum 
j- strings with the rapidities Za' ■ This state will have a total spin 

Sz = '^ rijVj - ^ Ukk > 0, (50) 
j k 

and energy 

^ = EE^'^(^')-EE^^-("''')' 4 e (-00; +00) \ A,, z^'^'gA,-. (51) 

k cx.=\ j a=\ 

The BAE poses restrictions on the possible combinations of excitation parameters {z^} and 
{zq'-' }. We miss here details of the treatment of this problem, which would involve the analysis 
of Eqs. (fTU)) in the next to leading order H15() . i.e., the 1/A^ corrections to (fT^. This is based 
on the known methods \\M\ and will be reported elsewhere. Instead we bring here some 
results for lowest-lying excitations with |Mj — Mf\ < 1. 

• Mi = Mf] The particle-hole excitations: 

E = -ei{z'''^)+ei{z^), = 0. (52) 



P\{x) > 0, p^ix) > 0, 
pI,{x)>0, pI,{x) = 0, 
P^j{x) = 0, j 7^1, mi, 
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• Ml = - I, Mi = Mf, i>l; The two-hole excitations: 

E = -ei(zf'') - = 1. (53) 

• Mjni = M^_^ — 1, Mi = Mf, i / mi; The particle-hole excitations: 

i? = -ei(z^'^)+ei(zi), S, = l. (54) 

• Mi = M^ - I, Mmi = M^^ +1, Mi = M^, i / 1, mi; The two-hole excitations: 

E = -ei{z'l'^)-ei{z^^^), = 0. (55) 

• Mi = M^ + 1, Mmi = M^^ - 1, Mi = M^, i / l,mi; The two-particle excitations: 

E = ei{zl) + eiizl), S, = 0. (56) 

All the rest excitations have |Mj — M-^\ > 1 for some i. Then, it is possible to have some 
combinations of excitations given above, as well as a new ones. For example, when po > 
— Mmi > ^ (> 1), Mi = M?, i / mi, together with other possible excitations it appears 
a new branch of two mi- hole excitation with the energy 

E = -Emi (^^1 '^) - Smi (2^2 '^)- (57) 

Notice that in any case the particle-hole excitation ()52() can be repeatedly excited: they can 
appear in any number as well as in combination with any other excitations. 

Besides these branches of the spectrum, there also exist bound states of magnons whose 
energies are given by ej, j = 1,2,.., mi — 1, (|45|) . These bound states are single-parametric, 
spin-singlet elementary excitations. Indeed, to construct a single bound state of s magnons, we 
must choose 

lk = Sk,s, pi{x) = 0, X ^ Ai, pi{x) = 0, X £ Ai, p^^=0. 
Then the j = 1 and j = rrii eqs. (|15|) acquire the form 

hix) = pi{x) + pi{x) + Til * Pi{x) + Timj * Pnn{x) + j;^Eis{x - Zs) (58) 

hmiix) = -pmiix) + Tmil * Pl{x) + Tminii * Pmi{x) + ^Tmis{x - Zg), 

where Zg parameterizes the excitation. From the last equation it follows that this configuration 
has a zero spin, 

Sz = Vl+ Vmi - sis = 0. 

It can be shown that H58() has a unique solution for Zg S (— oo; +oo). It constitutes elementary 
excitation with the energy 

£ = es{zs). (59) 

As seen from H'28() and (|47() . strings with j > mi have vanishing energies and become important 
only when one classifies excited states with respect to the total spin Sz. 
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Due to the properties of spectral functions e (|^^ - (|47() . excitation energies can take any 
positive values close to zero, then excitation parameters {z^^ and {z^^ come close to Fermi 
points or Zg ±00, i. e., the spectrum is massless, as it was expected. 

In the region 7r/2 < A < vr, we use transformation ()19|) which produces a mapping of this 
region onto the previous one. From (|19jl and due to the fact that Si{x) and di{x) in are 
even functions of x whereas hj{x) and eo{x) in (|^^ are odd functions of x, one gets that all of 
our solutions and results remain valid in this region also, when we change A by vr — A and the 
variable x by —x, in the corresponding functions. 

In the second case (j2| of imaginary A we deal with the spectral system (|4H) . It can be easily 
shown that 

TT 

/sinh[(n — 1)A1 x , „ , 

-— etWdij>0 for n>2 60 

cosh[(n — IjAJ — cos(v9 — ip) 

— TT 

obeys (|1T|) if £1(93) is solution to the integral equation 

TT 

ei{ip) = ^eoM + j r{ip- ^)eti^)d^, (61) 

— TT 

wher the real positive valued kernel r{ip) = -^Yl 2cosh(Afc) ^^^"^^^ periodic with 27r. As in the 

k 

previous case, from this integral equation we conclude that £1(97) takes negative values in some 
region A C (— cxd,0). The vacuum is formed by 1-strings filled with rapidities from Ai; the 
vacuum densities satisfy 

p?M = 0, ^0Ai, p? = 0, v^G Ai, 

Pniy^) = 0, ((C E [-vr,7r], n = 2,3,.... 
By the integration of (|34() for n = 1 one gets for the vacuum configuration 

Ml + ^Mi = N, Ml, Ml / 0, M„ = 0, n>2. (62) 

This means that our vacuum in this case has spontaneous spin 7^ 0. The band- filling is 
1^ < 1/2 (the exact value of band filling can be extracted from the exact solution of (|6H) and 
(HI) only ). 

The excitation spectrum includes massless 1-string particle-hole excitations with energies 
f = J]£i(Vp)-^£i(V',^), V'p^Ai, VgGA, (63) 

P Q 

and massive breathers of order n, n >2 with energies 

£n = Y.en{i^p). (64) 
p 
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4 Some Special Limits 



4.1 Ji = 0, J2^0 

This case is also in the region It corresponds to the hmit 9 — > itioo, when Ji = and 
J2 = =F^tanA. Consider the integral equation H48|) for 9^1 (we are changing 9 i9). By 
the definition (|22jl . eg depends on the combinations x ^9/\ and is essentially non-zero in the 
vicinity of ziz9/X. One can expect the solution si has the same property and expand ()48() near 
to these points: the mutual influence of these regions will be of order 0(e~'^' Then (|48j) goes 
to the following two decoupled integral equations: 

vrtan A/A 



ef(x) 



cosh(^x) 

vrtan A/ A 
cosh(^x) 



R{x - y)ef{y)dy, 



(65) 



ef >0 



+ 



R{x-y)e{{y)dy, 



£1 >o 



where e{{x) = ei{x - 9/\), ef (x) = £i{x + 9/X). 

By their forms, Eqs. (|65() are the same as corresponding integral equation for the spectral 
function £1 of Heisenberg XXZ model (see for example |26|): from this naive point of view, here 
we have two samples of XXZ chains with different signs of coupling constants. 

It is natural to try to confirm this double chain structure starting from the Hamiltonian, 
which on the other hand would be seen a self-consistency check for our results. In the case under 
consideration, the Hamiltonian takes the form 



N 

^=EE{(-i)V2 



j=l s=0,l 



Being rewritten as follows: 
1 ^ 



(66) 



-|- tan A i 



exp{(-l)^iAcj|^i„ J • ^j',s^j+i,s + ^-c- 



cos A 



(67) 



it reveals the underlying double chain structure, which is more evident under the following 
gauge-like unitary transformation: 



H 



1 



E 



cos A 



N>n>k>l 



T. O'l 1 



where 



H 



N 



XXZ 

s, open 



N>n>k>l 



N-1 

E ^M'^i+i,^ 
1=1 



+ + 7j COS A al,af_^^^,. 
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So the system can be described by the unitarily equivalent Hamiltonian 




H^^^ + (e ''V+ o^i;o + h.c.) - Hr - (e 





a+ + /i.e.) V (68) 



By the further apphcation of Bethe Ansatz technique, one can diagonahze this Hamiltonian. 
Then one parameterizes the eigenvalues of (|66p by means of two sets of complex rapidities {u*}, 
as attached to the chains H^^^, s = 0, 1. The eigenvalue problem leads to the Bethe Ansatz 
equations 



The phase factors reflect the twisted boundary terms in the Hamiltonian H. On the other 
hand, these phase factors reflect an interaction between the two sets; one could assume, that the 
s = spins are inside a "gauge field" generated by the s = 1 spins and vice versa. The energy 
eigenvalues are given by the formula 



which does not affect the mentioned interaction explicitly. 

Let us mention that the same equations one would get from BAE (jSJ. Indeed, in the limit 
9 ^ 1, {9 ^ i9), the roots of BAE (j^J can be divided into two groups gathered close to the 
points ±9/X: = u + 9/X and v} = u — 9/\. Taking the limit ^ cxd in the BAE one 
would get (inH). 

It is intuitively clear that, in the thermodynamic limit N ^ oo, boundary terms does not 
affect the bulk behavior and the Hamiltonian H66|) will have the same critical properties as two 
non-interacting XXZ- chains with usual periodic boundary conditions and opposite signs of 
coupling constants. This is exactly what we had in (|65|1 : the s = and s = 1 chains correspond 
to L and R sectors respectively. 

The interaction effect is of the next to leading order (actually of order l/N"^). It can be 
studied using the known finite size correction methods [301 1^ 02] . The idea is to consider two 
XXZ chains, each with additional phase shifts in the Bethe Eqs. (|in|l . The phase shift of one of 
the chains is described by the total spin of another chain. This requires rather long calculations. 
We bring here only the result for the energy difference between the ground state and the one 
with {n^) particles added at the right (left) Fermi points of the s = 0, 1 chains: 




(69) 




(70) 



2-K 
4/V2 




3 



1 



(71) 



where vq is the Fermi velocity and ^ is the dressed charge introduced in |^ 




(72) 
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and K^'^ are some integer quantum numbers. The interaction parameter 5 comes from the 
shifted phases. It is defined as 

Ni-s^s = {^^1--} ' « = 0' 1' (73) 

where {x} is the fractional part of x. Quantum numbers K are unimportant here, because are 
factored from the remaining interaction. We also used notations Ns = nf' + nf', Jg = — 



s 



Thus we get the "effective Luttinger liquid" description of our motel 

A similar integrable model containing two species of XXZ chains with the same signs of 
coupling constants is discussed by Schulz and Shastry in j^Hl • There the equations equivalent to 
H69() are derived in terms of momenta of Bethe excitations. The fact that the hopping parameters 
are of the same signs makes the physics different from our case. 

Following Ref. jl7j . we finish this section with the bosonization study of (|66|) . In order to 
obtain a right continuum theory, we first apply the Jordan- Wigner transformation along the 
zigzag links. The continuum fermionic Hamiltonian is obtained by the replacement Vn.s — *■ 
^/oq {Rs{x) + (— l)"'Ls(x)), X = naQ, where 'ipn,s is the discrete Fermion field, R{x) and L{x) are 
the continuous right and left Dirac fields respectively, oq is the lattice spacing. In the leading 
order we get: 



H = Hn + H., 



•inti 



Ho = vol [ dxJ2 {Rtix)dRsix) - Lt{x)dLsix)) , 

s 

Here = RfRs, Jl = LfLg are the right and left fermion currents respectively, vq is the 
Fermi velocity and g = 2cotAtanh0 is the coupling constant. Then, using the Bosonization 
prescriptions, we arrive at the bosonic Hamiltonian 

H = vojdx[{l- + (9$i)2) + (1 + |_)(n2 + {d^^f) + ^-{liod^i - HiScDo)} . (74) 

Here we introduced Bose fields ^s{x) and canonical conjugate momenta ns(a;) with the commu- 
tations 

[^s{x),Tis'{y)] = iSss'S{x - y), 

In deriving H74|) we have dropped strongly irrelevant interaction terms like cos(4-y/7r ^^(x)). It 
is not hard to check that in the first order in small A the Hamiltonian (|74)) gives the energy 
difference equivalent to (f7T|) . The latter allows also to extract an exact values for Fermi velocities. 
The canonical transformation 

_ 1 9$o + ni _ 1 9$i + no 



V2 ' ~ V2 
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brings ()74|) into the Gaussian form 

H = voK+ j dx {Pi + dQl) + vqK^ j dx + dQ'i) . (75) 

We thus see that the interaction between the two chains effects on the renormahzation of Fermi 
velocities only. 

4.2 Ji = J2 

This special limit is a crossing of two regions (pQ) and (jSJ- It takes place when A — > and — > 
with finite ratio ^ = const. In order to consider this limit, it is convenient to redefine 9 i\9 
in the region (0), and A — > zA, ^ — > —i\9 and if — > A(/?, in the region (here is the rapidity 
variable introduced in the section (2.2)). In both cases, we will have 

Ji = J2 = -i^. 

Consider this limit within (^. Here one has po ^ 00 and vi = mi 00. This makes the 
system of coupled spectral equations H25() simpler; the pseudo-energies ej become non-negative 
for j > 2, like in the case of H41() . Then (|24() for j = 1 goes into the integral equation for £i{x) 
as follows: 

ei{x) = AoKx) + j To{x-y)ei{y)dy, (76) 

ei<0 

where A = 2 (l + 1/9"^) , Tq{x) = ^ (^7^^, h{x) = - ■ ^lie rest for j > 2 

can be expressed in terms of the solution of (f7B)) through where 



00 



G,(x)= / ^e---e-(-'-i)l-l ^ ^^-^'"^^ 



^ ' 27r 27r j;2 + (j _ i)2 • 



-00 



Let us mention that one in principle could take the limit in H48() . Then of course the obtained 
integral equation will be equivalent to (|7H|) : instead of our case, it will contain an integral over 
the region ei > 0. 

The same equation could be found in the region ((2j. Indeed, the rescaling 93 Xip destroys 
the periodicity — > 93 + 27r and the new 93 varies from —00 to +00. Taking this in account, 
one can make sure that (|39j) for n = 1 goes into (|76|1 . This proves that we indeed get the same 
pictures coming from two different regions and ©• 

From (|76j) one can conclude that ei takes negative values in some region Aq C (— oo,0). 
Thus the vacuum in this case is filled by 1-strings with the rapidities from this region; the 
corresponding vacuum densities with the properties 

p^{x) =0, X Ao, ^0(2;) =0, X G Ao (77) 

obey the integral equation 

b{x) = Pv{x) + Pv{x) + j Tq{x- y)p^{y)dy, (78) 

Ao 
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where b{x) = ^ (^x+e)^+i ~^ (x-e)'^+i ■ This integral equation together with ((77|) uniquely 
defines the vacuum densities. The vacuum configuration have properties very similar to the one 
from the section (2.2), given in ()62() . 

Another phase transition occurs when one in addition sends 9 — > ±00. It corresponds to Ji = 
J2 = 0. This case can be treated similar to the special limit from the previous section. It can be 
checked that (|76|) splits into two integral equations which are equivalent to corresponding spectral 
equations of isotropic XXX-chain. The two integral equations will diff'er by opposite signs 
of the free term (disturbance function), which means that the corresponding isotropic chains 
are one ferromagnetic and another antiferromagnetic. The same follows from the Hamiltonian 
straightforwardly, when one puts Ji = J2 = 0. This is another consistency proof of our results. 



5 Summary and Conclusion 

To summarize we have applied the method introduced by Gaudin in and Takahashi and 
Suzuki in [2^] to the exactly solvable two-leg ladder model with zigzag like interaction con- 
structed in We have established two distinct zero temperature phases, corresponding to 
two regions in coupling constant space IJ1/J2I < 1, O and IJ1/J2I > 1, (O. We considered 
the model without magnetic field. When IJ1/J2I < 1, we found the model has gapless excita- 
tions above the antiferromagnetic ground state formed by filling all the possible 1-strings (even 
magnons) with rapidities from some region Ai C (—00, 0) and mi-strings (odd magnons) with 
rapidities from (—00, +00). So this phase is conformal. In the region @, we have seen the 
model possesses a ferromagnetic ground state with gapless particle-hole excitations and massive 
breathers of length n > 2. The intersection of these two phases |Ji| = IJ2I was considered in 
the last section. In the last section, we have considered also the special limit when Ji = 0. We 
have shown that in this case our spectral equations are equivalent to analogous equations for 
two XXZ chains with opposite signs of coupling constants. This double chain structure was 
confirmed by a unitary transformation of the Hamiltonian. We would like to mention that 

It naturally arises a question about the scaling dimensions of fields when the conformal phase 
occurs. These calculations could be done by use of the method introduced in |36j and [Slj . 

The problem which did not allow to get a precise values of Fermi points, Fermi velocities, 
momenta of states, as well as conformal dimensions is to solve a Fredholm type equations as 
(|15|) and (|HT|) are. As we already mentioned, the perturbative Wiener-Hopf method could of 
course be effective, though it does not usually give an explicit solutions. The numerical analysis 
is also appropriate. We plan to present this analysis and related results elsewhere. 
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